Anomalously slow phase transitions in self- gravitating systems 
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Kinetics of collapse and explosion transitions in microcanonical self-gravitating ensembles is analyzed. A 
system of point particles interacting via an attractive soft Coulomb potential and confined to a spherical container 
is considered. We observed that for 100-200 particles collapse takes 10 3 - 10 4 particle crossing times to 
complete, i. e., it is by 2-3 orders of magnitude slower than velocity relaxation. In addition, it is found that the 
collapse time decreases rapidly with the increase of the softcore radius. We found that such an anomalously long 
collapse time is caused by the slow energy exchange between a higher-temperature compact core and relatively 
cold diluted halo. The rate of energy exchange between the faster modes of the core particles and slower-moving 
particles of the halo is exponentially small in the ratio of the frequencies of these modes. As the softcore radius 
increases, and the typical core modes become slower, the ratio of core and halo frequencies decreases and the 
collapse accelerates. Implications to astrophysical systems and phase transition kinetics are discussed. 
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I. INTRODUCTION 

Many groups of stellar systems have highly universal struc- 
tures despite the apparent differences in their history and en- 
vironment Illy]. These universal structures are thought to 
have arisen as a result of relaxation towards equilibrium or to 
otherwise long-lived states. A comparison between the age 
of the stellar systems with universal features and correspond- 
ing collisional relaxation times reveals, however, that several 
types of stellar systems, such as elliptical galaxies, have not 
existed long enough to be collisionally relaxed |1]. Other 
than collisional types of relaxation, such as "violent relax- 
ation", or phase-space mixing caused by strong gravitational 
field fluctuations |3], have been suggested to explain this ap- 
parent contradiction between the timescales. Yet a full un- 
derstanding of the kinetics of relaxation in naturally occurring 
self-gravitating systems systems is still lacking. 

A number of fairly idealized models have been analyzed 
to understand the nature of equilibrium and transitory states 
of stellar systems. A well-studied example is an ensemble 
of self-gravitating particles particles with a sufficiently short- 
range small distance regularization confined in a container. 
That system exhibits a gravitational phase transition between 
a relatively uniform high-energy state and a low-ener gy s tate 
with a core-halo structure 0, H S S S S ED- EH 
fl^l . During such a transition in a microcanonical ensemble 
the system undergoes a discontinuous jump from a state that 
just ceases to be a local entropy maximum to a global entropy 
maximum state with the same energy but different tempera- 
ture. A transition from a high-energy uniform state to a lower- 
energy core-halo state is usually called collapse. The reverse 
transition during which the core disappears is often referred to 
as explosion. 

It has been recently observed in molecular dynamics (MD) 
simulations 1 15] that a typical timescale for such gravitational 
transitions is paradoxically large, for a system of 125 - 250 
particles being of the order of 10 3 relaxation times for a col- 
lapse and 10 2 relaxtion times for an explosion. The relaxation 



time, t r = FP^N 1 / 2 /hi N, is the timescale of typical par- 
ticle velocity thermalization which proceeds mostly via soft 
Coulomb collisions (l|. It was also observed in Ref. 111511 that 
the density relaxation, such as the formation of the core, ad- 
vances relatively fast, while the evolution of the kinetic energy 
or temperature proceeds noticeably slower 11511 . 

In this paper we undertake a more detailed study of kinetics 
of collapses in self-gravitating systems, in some way com- 
pleting the investigation initiated in [15]. The structure of the 
paper is following: After this introduction we briefly outline 
the simulation setup and present the results for collapse kinet- 
ics. A section analyzing a slow core-halo energy transfer as 
the bottleneck of system relaxation follows. Conclusions and 
discussion of the results completes the paper. 



II. SIMULATIONS 

We consider systems consisting of N = 125 — 250 identi- 
cal particles of unit mass confined in a spherical container of 
radius i? with reflecting walls. The Hamiltonian of the system 
reads as 
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where ro is the soft core radius. Along with the physical 
units, we use the standard rescaled units (as discussed, e. g., 
in Ref. [4]). For energy e, temperature 9, distance x, and time 
r they read as 
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FIG. 1: Plots of entropy s(e) (solid line) and temperature r(e) = 
de/ds (dashed line) vs energy e for a system with a gravitational 
phase transition and a softcore radius xq = 0.005. 
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FIG. 2: Plots of the relative values of (from top to the bottom) the 
number of core particles N' c (t) (blue), core temperature c (t) (al- 
most coinciding with N' c (t), but with stronger fluctuations, dotted 
black), and total temperature Q'(t) (red) vs r for a collapse in sys- 
tem with e = -0.5, N = 125, and x = 0.005 



Expressed in these rescaled units, equilibrium properties of 
self-gravitating systems become universal. The velocity re- 
laxation, assuming that it is caused mostly by the soft col- 
lisions, is expected to be universal in terms of time r ve \ — 
t In N/N where the factor N / In N is proportional to the 
number of crossings a particle needs to change its velocity by 
a factor of two. 

The phase diagram of the system is presented in Fig.^ see 
also Refs. 181 fl21 [1511 . High- and low-energy branches termi- 
nating at the energies e co u and e exp i correspond to the uniform 
and core-halo states. The collapse and explosion energies are 
e co ii ~ -0.339 and e exp i « 0.267 for x = 5 x 10~ 3 . 

Each MD run was initiated with a configuration in which 
the particles were seeded according to the corresponding e co u 
or e exp i equilibrium (metastable) density profiles, and the ve- 
locities were assigned according to the Maxwell distribution. 
A more detailed description of the MD simulation procedure 
is presented in 11511 . 

To reveal all facets of gravitational phase transitions in the 
most informative way, we consider the following parameters: 

• Temperature O as the indicator of the advancement of 
a phase transition as a whole. When 6 reaches the tar- 
get phase equilibrium value, all other system parame- 
ters come to equilibrium as well and the phase transition 
is complete. 

• Number of core particles N c as the measure of density 
relaxation. 

• Temperature of the core 9 C which is proportional to the 
average kinetic energy of the core particles. Deviations 
of 9 C from 6 quantify the temperature gradients occur- 
ring during a phase transition. 

We do not list here any parameters which characterize 
the velocity relaxation: As it follows from the definition of 



T V ei OL an d as observed in simulations [ 15], the velocity dis- 
tributions in both core and halo become thermalized within 

Tvel ~ 1. 

Results, averaged over 4 runs, for the temperature, core 
temperature and the number of core particles for a collapsing 
system are presented in Fig. [2] Time evolution of these param- 
eters is plotted in terms of relative variables O'(r), Q' c (t), and 
N' c (t) which are denned as 9'(r) = [9(r) - 9(u)]/[9(c - 
h) — 0(tt)]. The values 9 (it) and 9(c — h) correspond to the 
uniform and core-halo states in equilibrium. 

As in Ref. 1 15], we observe that a collapse in a system with 
N = 125 — 250 particles and x = 0.005 takes about 10 3 
velocity relaxation times to complete. It also follows from 
Fig. |2] that the growth of the core is significantly faster than 
the relaxation of the average kinetic energy: The core reaches 
half of its equilibrium size only in about 5 velocity relaxation 
times, while temperature relaxes to half-way only in 110 ve- 
locity relaxation times. In addition, we can conclude that the 
core temperature is evolving almost synchronously with the 
number of core particles, i. e., noticeably faster than the total 
temperature of the system. 

The results for the collapse in an otherwise identical system 
but with softcore radius twice larger, x$ = 0.01, are presented 
in Fig. [3] It follows from the comparison between Figs.|2]and 
[3]that while the initial stages of relaxation are not affected by 
the change of the short-range potential, the overall collapse 
proceeds much faster for larger xq. 

The above numerical results suggest that after a rapid ini- 
tial core growth, which takes few velocity relaxation times, a 
further evolution is hindered by some slow process that is es- 
sential for the completion of the phase transition. Kinetics of 
this slow process strongly depends on the short-range part of 
the interparticle potential. It follows from Figs. [2] and [5] that 
a collapse comes to its completion only when the tempera- 
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FIG. 3: Same as in Fig.|51but for x = 0.01 



tures of the core and the halo become equal. Hence, it seems 
natural to assume that the bottleneck process is the energy, or 
heat, exchange between the core and the rest of the system. 
The rate of this heat exchange depends on the structure of the 
core which in turn is determined by the potential softening. In 
the next section we consider the heat exchange between the 
core and halo in more detail. 



III. CORE - HALO ENERGY EXCHANGE 

To analyze the energy exchange between the core and halo, 
let us first examine the motion of a core particle. In a system of 
reference of the center of mass (CM) of the core we expand the 
potential energy terms in the Hamiltonian, Eq. Q, in powers 
of Vij I tq and arrive to the harmonic oscillator Hamiltonian, 
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Hence the motion of core particles relative to the core CM is 
characterized by harmonic oscillations with the frequency 



(4) 



It is interesting to note that the frequency of the motion of 
a particle in a uniform self-gravitating sphere of N c particles 
and radius ro with the bare gravitational (— l/r) interaction is 
also given by Eq. Q. Since the core radius is roughly equal to 
fo 11511 ■ both the bare interaction and the "very soft" potential 
frequencies are essentially the same. 

While the motion of the core particles is bound by grav- 
ity, the higher energy halo particles can be viewed as free and 
being confined by the container walls only. Hence, the inter- 
action between a core and a halo particle can be approximated 
as an interaction between a point mass on a rectilinear trajec- 
tory and a harmonic oscillator. The first relevant term in the 
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FIG. 4: Sketch of a core - halo particle scattering event. 



multipole expansion of this interaction is the monopole-dipole 
term 
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Here Xi and Zj are the coordinates of core particles, v and p 
are the velocity and the impact parameter of the halo particle, 
respectively (Fig. 0). The monopole-monopole term, corre- 
sponding to the interaction between the halo particle and the 
core CM, is irrelevant to the internal motion of the core parti- 
cles. 

With the introduced approximations the core-halo energy 
exchange becomes physically identical to the well-studied 
process of energy exchange between the vibrational degrees 
of freedom of bound states and fast free particles in a plasma. 
Naturally, such molecular processes are usually considered 
in quantum mechanical terms. Following a standard text- 
book 1 17] we start with the first-order term of the interaction 
representation expansion of the perturbation of the oscillator 
wave function during a complete single collision, 
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and arrive to the following expression for the probability |C| 2 
of the excitation from (or de-excitation to) the ground state of 
a harmonic oscillator, 
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Here Kj are the McDonald functions (Modified Bessel func- 
tions of the Second Kind), the factor N 2 appears since the 
interaction potential in Eq. (0 is the sum of N c identical 
terms and the probability is quadratic in V(t). Multiplying 
Eq. |7} by the transferred energy %u) c , we get rid of the quan- 
tum constants and obtain an expression for the typical energy 
exchange during a collision between a core and a halo particle, 
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The same expression can be obtained by completely classical 
analysis considering the energy transfer during forced oscilla- 
tion 11811 . The last limit in Eq. ^ is taken since p/v, which is 
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of order of a halo particle crossing time R/v, is much larger 
than a period of oscillation of the core particle l/ui c . It fol- 
lows from Eq. l|8) that the rate of energy exchange between 
the core and halo particles is exponentially small in the ratio 
of typical frequencies of their motion. 

To obtain the rate of energy transfer between the core and 
all halo particles per unit time, one needs to average Eq. (|8} 
over impact parameters and velocities, 



2irpdp / vW M (v)nSE(p,v)d 3 v. (9) 
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Here the lower limit of integration for the impact parameter is 
set equal to the core radius for the dipole approximation to be 
correct, Wm{v) is the Maxwell distribution, n = 3N/ (4irR 3 ) 
is the particle density. Using that tq -C R and evaluating 
the velocity integral in the steepest descent approximation, the 
following expression is obtained 
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Assuming that the typical time of core - halo relaxation t c -h 
is roughly equal to the ratio of the total transferred energy 
to the rate of the transfer and using Eq. @, we arrive to the 
following expression for t c ~h, 



tr.-h 



9 2 /3 (a;o / nc )l/6 



7r\/3T' 



■ exp 



1/3' 



(11) 

Here Ae is the total transferred energy in rescaled units and 

n c = N c /N. 

To evaluate the numerical value of t c -h we consider the ex- 
ample from Fig.|2]with the following parameters: e = —0.5, 
x = 0.005, n c » 0.2, 9 « (e^ - « 1.2 (which 

is the average between the initial and final halo tempera- 
tures) lll5ll . We assume that the total change in kinetic en- 
ergy of the system during the collapse is a result of the energy 



transfer from the core and set Ae 
these parameters we obtain 
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for the relaxation time. Given the number of approximations 
used in obtaining Eq. 1^3, the agreement with the simula- 
tional result for the complete phase transition time t^V ps 
27000 (see Fig. [2) is surprisingly good. The agreement is even 
better for systems with the same energy e = —0.5 but larger 
softcore radius, xq = 0.01, where n c ss 0.22, 8 ps 0.83, 
and Ae ps 1.5. The theoretical estimate, Eq. il It . yields 
r c -h ~ 3600 which is only very little below the simulational 
result re 3800. 

The most significant contributions to underestimating the 
value of Tc-h are the following: 



• All collision were considered complete, i. e., the inte- 
gral in Eq. l|6} had infinite limits. This is certainly not 
true for a confined system especially for large impact 
parameters. 



• The energy exchange between the core oscillations and 
halo particles was always considered in one direction, 
i. e., from core to halo. In reality, this is only true for 
slow halo particles. Close to equilibrium, the exchange 
becomes progressively mutual with increasing number 
of fast halo particles losing their energy to the core vi- 
brations. 

Other factors not taken into account in our estimate, 
but possibly affecting the core-halo energy relaxation are: 
Quadrupole and higher-order terms in the potential expan- 
sion, overlapping collision, more complex dynamics than os- 
cillation and rectilinear motion of the core and halo particles, 
higher than one photon processes, or higher order terms in 
the perturbation expansion to mention some. Yet we believe 
that our relatively simple approach presented above captures 
the essence of the core-halo energy relaxation and has semi- 
quantitative predictive power. 



IV. DISCUSSION AND CONCLUSION 

In the previous two sections we obtained the following re- 
sults for the kinetics of collapse from the uniform to the core- 
halo state in self-gravitating systems: 

1. The molecular dynamics simulations revealed that in 
a system of 100-200 particles the collapse time is by 
two-three order of magnitude longer than the velocity 
relaxation and strongly depends on the short-range part 
of the interaction potential. 

2. It was found that the non-equilibrium feature with the 
slowest relaxation time is the temperature difference 
between the core and the halo. In contrast, such pa- 
rameters as number of core particles and core kinetic 
energy evolve relatively fast. 

3. A mechanism similar to the vibrational-translational 
relaxation in plasmas was suggested for the core - halo 
energy exchange. For this mechanism we show that the 
core-halo thermalization is exponentially slow in the 
ratio of typical frequencies of the motion of core and 
halo particles. Despite several rather strong approxi- 
mations used in our analysis, a theoretical estimate for 
the relaxation time is in a good agreement (less than by 
a factor of 2 off) with the simulation results. 

So far nothing has been said about the reverse to collapse 
transition, i. e., explosion, illustrated in Fig. [5] It follows from 
this Figure that similarly to the collapse, the kinetic energy is 
the slowest-evolving quantity while the number of core par- 
ticles and core kinetic energy lead the explosion. Since the 
fastest particles leave the core first, the core is always colder 
than the rest of the system. At the end of explosion the core 
becomes just a single cold particle indistinguishable from any 
other particle of the system. This explains the sudden jump 
of the core temperature at the very end of explosion. Us- 
ing Eq. it i s straightforward to explain why an explosion 
is faster than a collapse: Since the explosion starts at higher 
energy than collapse, the corresponding 'explosion core-halo 
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FIG. 5: Plots of the relative values of (from top to the bottom) num- 
ber of core particles N' c (t) (blue), total temperature 0'(t) (red), and 
core core temperature O c (t) (black) vs r for the explosion in system 
with e = 0.5, TV = 125, and x = 0.005 



state' has noticeably fewer core particles (twice for the con- 
sidered here case) than the final collapse state. This reduces 
the exponential term in Eq. (^3. In addition, the total amount 
of energy which needs to be transferred between core and halo 
is smaller in case of explosion than in the case of collapse. It 
follows from Fig. ^ where temperature jumps at collapse and 
explosion points are proportional to the energy exchanged be- 
tween core and halo. 

Are these results applicable to self-gravitating systems with 
other types of short-range regularization? For systems with 
continuous potentials at low-energies the particle motion near 
the equilibrium position are harmonic oscillations with a fre- 
quency roughly given by Eq. @ where ro is of order of 
the core radius. Examples include ensembles with Fourier- 
truncated Coulomb potentials ffl or exchange interaction in 
systems with phase-space exclusion |6]. For potentials with 
singular short-range part such as hard-core repulsion, the mo- 
tion of core particles is discontinuous, yet a typical inverse 
timescale of such motion is given by Eq. as well. This 
follows from the fact that the expression Eq. (|4j can be ob- 
tained by dividing a typical core particle velocity y/N c /ro by 
the core radius ro. The analysis and conclusion made in the 
Section III are based on a general principle that for a perturba- 



tion of a fast system by a slow one, the rate of energy transfer 
between these system is exponentially small in the ratio of 
their frequencies. And since the motion of core particles is 
always faster than the motion of the halo ones, we conclude 
that the results obtained here are qualitatively applicable to all 
collapsing self-gravitating systems independent of the form of 
the short-range cutoff. 

An important implication of the demonstrated here essen- 
tial non-isothermicity of collapsing system is about the ap- 
plicability of the Smoluchowski equation description to these 
systems Hill . Compared to the more complete Boltzmann 
equation with Landau collision term (often called in case of 
self-gravitating systems Fokker-Plank-Vlasov equation 1 16]), 
Smoluchowski equation offers a significant simplification: In 
the physically relevant case of spherically symmetric systems 
the solution depends only on one radial coordinate r. To make 
such description more realistic, one needs to take into account 
the non-uniform temperature field. This can be achieved by 
coupling the Smoluchowski and heat conduction equations. 
The non-isothermicity of evolution also indicates that the col- 
lapse in canonical ensemble must be radically different from 
its microcanonical counterpart. 

Finally few words about the astrophysical relevance of the 
obtained results. The structure and mere existence of the 
equilibrium core is a consequence of short-range regulariza- 
tion of the gravitational potential and the confining container 
which are the artifacts of the model. Hence the true equilib- 
rium core-halo states never occurs in stellar systems. How- 
ever, the observation that the collapse progress in hindered by 
the core - halo energy exchange which is exponentially slow 
in the ratio of the typical frequencies of motion of the core 
and halo particles remains applicable. Hence it is possible 
to interpret the states of the systems such as globular clusters 
whose age is noticeably larger than corresponding velocity re- 
laxation (and therefore, initial core formation) time 1 1 ] as the 
transitory long-living states similar to core-halo thermaliza- 
tion states observed in this study. 
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